Growth rate distribution of NH 4 C1 dendrite and its scaling structure 
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Scaling structure of the growth rate distribution on the interface of a dendritic pattern is inves- 
tigated. The distribution is evaluated for an NH4CI quasi-two-dimensional crystal by numerically 
solving the Laplace equation with the boundary condition taking account of the surface tension 
effect. It is found that the distribution has multifractality and the surface tension effect is almost 
ineffective in the unscreened large growth region. The values of the minimum singular exponent 
and the fractal dimension are smaller than those for the diffusion-limited aggregation pattern. The 
Makarov's theorem, the information dimension equals one, and the Turkevich-Scher conjecture be- 
tween the fractal dimension and the minimum singularity exponent hold. 

PACS numbers: 05.45. Df, 89.75.Kd, 68.70.+W 



I. INTRODUCTION 

The dendritic pattern is observed in various systems, 
such as crystallization[l| and electrodepositionj^, H|. It 
has been one of the most typical and ubiquitous in non- 
linear and nonequilibrium physics. A dendritic pattern 
has a stem with the tip growing stably and steadily, with- 
out splitting. Countless sidebranches grow behind the 
tip due to noise [4j and instability of a flat interface @. 
These sidebranches compete with the ones around them. 
A longer sidebranch screens off growth of the shorter 
ones around it. The competition repeats on various 
length scale, complicated and hierarchical structures are 
formed [(| and the pattern becomes fractal as a whole. 

The growth process of a dendritic pattern is mainly 
dominated by diffusion and anisotropy. The dendritic 
pattern is often compared with the diffusion-limited ag- 
gregation (DLA)Q pattern, for which the growth pro- 
cess is also dominated by diffusion. However, there is 
no anisotropy in the DLA growth process and the DLA 
pattern is formed through repeat of tip-splitting. It has 
been reported for an electrodeposition experiment that 
a transition between dendritic pattern and DLA pattern 
is observed as the electrolyte concentration and applied 
voltage are varied @, Furthermore, the dendritic pat- 
tern can be formed artificially by introducing anisotropy 
into isotropic viscous fingering [8, 9] and, on the other 
hand, the DLA pattern by removing the anisotropy from 
an anisotropic crystal growth process |10j. 

The scaling structure and fractality of dendritic pat- 
tern have been interesting and important issues for quan- 
titative characterization. It has been reported that the 
area fractal dimension Df of a dendritic pattern with 
fourfold symmetry is 1.5-1.6 for a noise-reduced DLA 
simulation on a square lattice ll| and NET^Br crystal 
growth|12]. It is clearly less than that of the DLA, 
Df ~1.71. This is attributed to the fact that the tip 
is stabilized against splitting by anisotropy. Mathemat- 
ically the lower bound of the fractal dimension for the 
DLA on a square lattice is proved to be 3/2 [l3| ■ 

To characterize a pattern in more detail, the fractal 
dimension alone is insufficient. For the DLA, the growth 
rate distribution on the interface is found to have mul- 
tifractality in the cases of a simulation [13]. a crystalliza- 



tion experiment [l5| , and a conformal mapping modelflfjl- 
Il8j . In the present article, we investigate the scaling 
structure of the growth rate distribution, especially com- 
pared with that of the DLA. We evaluate the growth rate 
distribution numerically, and implement the similar mul- 
tifractal analysis for a dendritic pattern formed in NH4CI 
crystallization experiment, where the surface tension is 
effective within a certain length scale. 



II. EXPERIMENT 



We use an NH4CI solution growth dendritic crystal 
with well-developed sidebranches. The details of our ex- 
periment are described in Ref.[6J: An NH4CI aqueous 
solution saturated at approximately 40° C is sealed, with 
a nucleation seed left in it, in a Hele-Shaw cell, a narrow 
space between two glass plates. The thickness of the cell 
is 100/im. Then when the temperature is lowered (to 
approximately 30°C), the solution becomes supersatu- 
rated, nucleation takes place from the seed and a grow- 
ing crystal is observed. The direction of the tip growth 
is (100) in the supersaturated solution and the growing 
dendritic crystal has fourfold symmetry. Sidebranches 
grow behind the tip and perpendicularly to the stem, 
with small sub-sidebranches perpendicular to them. The 
image of the crystal is obtained by a microscope and a 
charge-coupled device (CCD) camera, whose resolution 
is 640 x 480 pixels. The image is binarized by an image 
processing. The binarized image of a crystal interface is 
shown in FigfT] 

The tip growth velocity Utip is 40 — 49^m/sec. In 
the case, sidebranches are well-developed within the 
shooting window, as shown in Fig[TJ However since 
the spacing between them is smaller than the diffusion 
length, they are still competing each other, not grow- 
ing independently. The diffusion length near the tip 
Zd = 2-D/utip, where D is the diffusion constant of NH4CI 
(2.6 x 10 3 ^m 2 /sec [Hj]) is > 100/im, and that near a 
sidebranch is longer than this. Therefore the growth is 
regarded as quasi-two-dimensional. 
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FIG. 1. Image of a dendritic crystal. The resolution is 
5.5/im/pixel. Its fractal dimension is 1.54. The region in- 
side the broken line is the cutoff region where the surface 
tension effect is neglected in our analysis. 



III. GROWTH RATE EVALUATION 

In principle it is a faithful method to the original data 
to evaluate the growth rate from the growth site area be- 
tween two successive images, as implemented for a DLA 
pattern[ll|. However, for a dendritic pattern, it is quite 
difficult to implement with satisfactory precision due to 
the limitation of the resolution and since the difference of 
the growth rates between the fast region (around the tips 
of the stem and longer sidebranches) and the slow region 
(deep inside the forest of sidebranches) is much larger 
than that for the DLA. Therefore instead, we evaluate 
the growth rate Pgr(r; n t) at a point rj nt on the interface 
by the gradient of the concentration field 4>(r), 

Pgr(r in t) ~ |V0(r)|, (1) 

where </>(r) is assumed to satisfy the Laplace equation, 
V 2 </>(r) = 0, outside of the pattern. This assumption is 
valid if the diffusion length is larger than the characteris- 
tic length scale of the system, for example, the tip radius 
of the stem or the average spacing of the sidebranch gen- 
eration. In this case the characteristic length scale is of 
the order of 1/xm, thus the condition is satisfied. 

The Laplace equation for the concentration field is nu- 
merically solved by the relaxation method on a square 
lattice, whose spacing is set to be the pixel size. The 
concentration is supersaturated far away from the inter- 
face (the saturated concentration at 40°C, 46g per lOOg 
water) and the Gibbs-Thomson boundary condition is 
imposed at the interface [l9l [20j . 

<P(r ia t) = M^ + d(BHr int )), (2) 
d{6) = do(l -cos[40]), (3) 



with stiffness and fourfold symmetry taken into account, 
d = 2.24 A[24j the capillary length, and K(r; nt ) the 
local curvature at r int , respectively. The growth angle 
9 is defined as the angle between the growth direction 
at ri n t and that of the stem. The curvature and growth 
angle are calculated by spline interpolation for the pixel 
data of the crystal interface. 

The range in which the growth rate takes value is vast. 
In the whole interface, the ratio of the largest growth 
rate to the smallest is more than 10 10 . Even within the 
region around the tip of the stem, the ratio is more than 
10 4 . 

When the surface tension effect is taken into account, 
there may be a lattice point where the curvature ra- 
dius is smaller than the lattice spacing. The curvature 
may vary abruptly from point to point around such a 
lattice point. If this situation occurs deep inside the for- 
est of well-developed sidebranches, around their roots, 
a large growth rate, of the same order as that of the 
tips of the stem or longer sidebranches, may be gener- 
ated. This is unnatural since the growth in the region 
is strongly suppressed. In order to appropriately take 
account of the surface tension effect and phenomenolog- 
ically circumvent the above unnatural situation, we con- 
sider three cases below: (i) The surface tension is com- 
pletely neglected. In other words, at the interface we set 
0(i"int) = 0o uniformly. Since the typical length scale of 
the system is of the same order as the lattice spacing and 
can be regarded as the length scale within which the sur- 
face tension is effective, this setting is reasonable. This 
case is labelled "Laplace". In this setting the growth 
rate distribution is the harmonic measure, (ii) A cutoff 
k c is introduced. If \n\ > k c , \k\ is replaced with k c . 
Here we set k c = 0.01 times the reciprocal of the lat- 
tice spacing. This case is labelled "cutoff(l)". (iii) The 
surface tension effect is taken into account only around 
the tips of the stem and sidebranches and is neglected 
deep in the sidebranches, in the region shown in FigJT] 
This case is labelled "cutoff(2)". The dependence of the 
results on how to choose k c for the cutoff(l) case and 
the cutoff region for the cutoff(2) case is very weak as 
long as k~ 1 is large enough and the cutoff region is wide 
enough, respectively, to suppress the generation of the 
unnaturally large growth rate. 



IV. MULTIFRACTAL ANALYSIS 

Let us consider that the interface of the dendritic crys- 
tal is covered by disjoint boxes of size e and let Pj(e) be 
the growth rate in the j-th box, 

Pi( £ ) = H Pgr (fint)- (4) 

rintGj-th box 

The rate pj (e) is normalized to be a probability measure, 

y~]j Pj(c) — 1, where N(e) is the number of boxes nec- 
essary to cover the interface completely. The generalized 
dimension D(q) is defined as|2l| 



where <fio denotes the saturation concentration(at 30° C, 
41g per lOOg water), d(9) the surface tension coefficient 



(5) 
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for g ^ 1, where Z(e,q) is the partition function 

N(e) 
3 



(6) 



For g = 1, 



D(l) = lim 

e->0 



Ef (e) Pj( £ ) lQ gPj( £ ) 

loge 



(7) 



Practically -D(g) is evaluated from the slope of log Z(e, q) 
for log e by least squares method. 

The singularity exponent a = a(q) and its fractal di- 
mension /(a) = f(a(q)) are obtained as functions of q 
by the Legendre transformation of the generalized di- 
mension -D(g)[12|: 



a{q) 



dq 



[(8 



/(a(g)) = ga(g)-(g-l)D(g). 



(8) 
(9) 



However, it is not practical to evaluate a and f{a) from 
© and ©, since it may produce relatively large numer- 
ical errors. Therefore instead, we adopt a direct evalua- 
tion method presented in Ref.[2^ described below. 

First let us construct a new probability measure 
fij(e,q) with parameter q frompj(e) as 



Mi(e>8) 



{Pj(e)} q 



Ef ej te( £ )> 9 

Then let us define £(e, g) and £(e, g) as 

C(e,8) = X^( e ' 9 ) lo g[ft( e )L 

i 

£( e >8) =^^'( e '8) 1 °s[/ i i( e '8)]- 



(10) 

(11) 
(12) 



From them a and /(a) are, as functions of g, given as 

CM) 



a(g) = lim 



/(g) = lim 



o loge 



e^O log e 



(13) 
(14) 



Practically they are evaluated from the slopes of £(e, q) 
and £(e,g) for loge, respectively, by least squares 
method. It is easy to show that the definitions ([13")) and 
(TTJ} satisfy the relations ((SJ and © by direct calcula- 
tion. 



V. RESULTS AND DISCUSSION 

We are interested in the larger growth rate regime, 
q > 0. For the pattern of FigfTJ the results of D(q) for 
the three cases are shown in Fig[5J The log-log plots 
of Z(q, e) against the box size e and least squares fit- 
ting for several values of q are shown in FigJ3] The box 
size is chosen from 4 to 80 pixels, taking the thickness 
of branches into account. There is a good agreement 
between the results of the three cases. 
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FIG. 2. Generalized dimension D(q), q > 0, for the three 
cases for the pattern of FIG[T] The increment Aq is 0.2. The 
error bars are obtained over 13 samples for the "Laplace" 
case. 
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FIG. 3. Log-log plots of Z(q, e) for the pattern of FIGQ] 
against the box size e for q=0,l,2, and 3. Note that the 
slopes mean (q — l)D(q). For visibility, the measure is not 
normalized to be a probability (Z(q = l,e) is the total sum 
of the growth rate). 



The results of the multifractal f-a spectrum for the 
three cases for the pattern of Fig [T] in the small a region 
corresponding to q > are shown in FigJH These spec- 
tra are evaluated by Eas.(|10 p -(|14 p and the plots of £(g, e) 
and £(gs e ) against loge are shown in FIGHJ There is a 
good agreement between the results of the three cases, 
except for £(g = 0, e). This disagreement is attributed 
to the surface tension effect, especially the contribution 
of the growth at points deep inside the forest of side- 
branches, with curvature |k| ~ n c . The relatively small 
difference of £(g = 0, e) between for the "Laplace" and 
the "Cutoff(2)" case and the agreement of £(g, e), q=l,2, 
and 3, for the three cases indicate that the surface ten- 
sion effect is almost ineffective in the unscreened large 
growth region. The fact that D(q) and /(a) take contin- 
uous values and depend on g or a means that the growth 
rate distribution has multifractality. 
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FIG. 4. Results of the multifractal f-a spectrum for the three 
cases for the pattern of FIG[T] The contact point with the 
line f(ot) = a gives the information dimension D(l) = a(q = 
1) = f(a(q = 1))- The increment Aq is 0.1 for q < 1 and 0.2 
for q > 1. The error bars, corresponding to q =0.5, 1.5 and 
3, are obtained over 13 samples for the "Laplace" case. 
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FIG. 5. Fitting of a and / for q =0,1,2, and 3. (a)plots of 
C(q,t) against loge. (b)plots of f(g, e). The fitting lines are 
shown for the "Laplace" case. 



Laplace 1.55±0.04 
Cutoff(l) 1.55±0.04 
Cutoff(2) 1.554=0.04 
DLA 1.713±0.003 



1.53±0.02 1.02±0.03 0.56±0.04 
1.53±0.02 1.03±0.03 0.56±0.04 
1.534=0.02 1.07L0.03 0.58±0.05 
~1.71 - 0.6654=0.004 



TABLE I. List of characteristic scaling exponents. The av- 
erages and errors are obtained over 13 samples. Note that 
Df is the fractal dimension of the area of the pattern and 
D(0) is the fractal dimension of the interface, on which 
the growth rate measure is defined, and that by definition, 
ttmin = D(q —> oo). The values for DLA are cited from the 
results of a conformal mapping modelfl7l. Il8(|. 



Some characteristic values for the scaling exponents 
are summarized in Table fflalong with those for DLA 
conformal mapping model [171 . Il8| for comparison. Both 
the fractal dimensions of area D f and perimeter length 
-D(O) are about 1.5, manifestly smaller than those for 
the DLA. This result agrees with the results of on- 
lattice simulation [11] and mathematics [13[. The small- 
est singularity exponent is obtained at the tip of the 
stem where the growth is most active and it is also 
smaller than that for the DLA. Furthermore, this agrees 
with the Turkevich-Scher scaling conjecture Df = 1 + 
amin[H,[13]) which argues that the fractal dimension de- 
pends only on the scaling behavior of the growth of the 
most active domain. Note that it is clear that the tip of 
the stem is the most actively growing domain for a den- 
dritic pattern, while it is reported that for the DLA, the 
most active growth domain is not the outermost tip(T8|. 
The information dimension D(l) is regarded as the frac- 
tal dimension of the active zone where the growth is not 
screened (28[. It is proved by Makarov(29j that exactly 
D(l) = 1 for the harmonic measure. Our results agree 
with the theorem, remarkably even though the surface 
tension is taken into account. 



VI. CONCLUSION 

We evaluated the growth rate distribution of an NH4CI 
dendritic crystal interface by numerically solving the 
Laplace equation and investigated its scaling property. 
The effect of the surface tension is taken into account 
as the Gibbs- Thomson boundary condition with some 
types of cutoff introduced based on phenomenologically 
plausible assumptions. We found that in the unscreened 
large growth rate regime the distribution has multifrac- 
tality and the surface tension effect is not essential. The 
fractal dimension and the value of the smallest singular 
exponent are smaller than that of the DLA and consis- 
tent with the previous results given in theory and sim- 
ulation. Our results agree with the Makarov's theorem 
for the harmonic measure, D(l) = 1, and the Turkevich- 
Scher scaling conjecture, Df = 1 + a m i n in spite of the 
surface tension effect. 
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